Submitted for publication in ApJ Letters 

Preprint typeset using WT^i. style cmulatcapj v. 12/16/11 



(N 

o 
o 



WHY BARYONS MATTER: THE KINEMATICS OF DWARF SPHEROIDAL SATELLITES 

Alyson M. Brooks 1 , Adi Zolotov 2 

(Dated: July 12, 2012) 

Submitted for publication in ApJ Letters 

ABSTRACT 

We use some of the highest resolution cosmological simulations ever produced of Milky Way-mass 
galaxies that include both baryons and dark matter to show that baryonic physics (energetic feedback 
from supernovae and subsequent tidal stripping) significantly reduces the dark matter mass in the 
central regions of luminous satellite galaxies. The reduced central masses of the simulated satellites 
reproduce the observed internal dynamics of Milky Way and M31 satellites as a function of luminosity. 
Including baryonic physics in Cold Dark Matter models naturally explains the observed low dark 
matter densities in the Milky Way's dwarf spheroidal population. Our simulations therefore resolve 
the tension between kinematics predicted in Cold Dark Mater theory and observations of satellites, 
without invoking alternative forms of dark matter. 
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1. INTRODUCTION 

There are fewer small satellite galaxies orbiting our 
Milky Way galaxy than predicted by the favored Cold 
Dark Matter (CDM) co smological model ([Moore et al.l 
119991: iKlypin et al.l fl999h . Theories often reconcile the 
discrepancy between the number of observed satellites 
and CDM predictions by invoking the suppression of star 
formation in low mass ga laxies, for example by UV heat- 
ing at reionization (e.g., lOkamoto et al.l [2008T ) . If only 
the most massive satellites form stars, this can bring the 
predicted number of luminous satellites down from thou- 
sands to tens, in line with observations. Even then a 
serious problem remains, as the most massive satellites 
predicted by CDM models a re still much too dense com- 
pared to what we observe (iBovlan-Kolchin et alJl201lL 



l20HlWolf fe BullocU2012HHavashi fe Chibal2012ft . The 
tension between the observations and the predictions 
of the CDM model have led some researchers to pro- 
pose exotic or alternative forms of dark matter to re- 
duce the central masses of satellites (jVogelsberger et all 
I2012t lLovell et al.l 12012b . However, the highest resolu- 
tion simulations available to date to study the inter- 
nal properties of satellites include only the dark mat- 
ter (DM) com ponent of galaxies, neglecting the effe cts of 
baryons (e.g.. IPiemand et al.l l2007t iSpringel et "all 12003 
IBovlan-Kolchin et alJl2012f >. 

Studies using high resolution simulations that in- 
clude baryons have shown that star formation and 
subsequent supernova feedback can a lter the DM 
struc tur e at the center of ga laxi es dNavarro et"al 
19961: iRead fe Gilmord 120051: iMashchenko et al 
2006|: iGovernato et al.l 120101 : iCloet-Osselaer et al 



2012 : PMaccio et al.l 120121:' 



2013 IGovernato et all l2012t 



iPontzen fc Governatc 
Zolotov et al.l 120121: 



Tevssier et all I2012D IGovernato et al.l (|2010|) and 



Pontzen fc Governatol (|2012f l demonstrated that high 



resolution alone is not sufficient to capture the effects 
of baryons on DM. Rather, it is necessary to adopt a 
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realistic treatment of star formation that restricts stars 
to form only in high density regions comparable to 
giant molecular clouds. The subsequent supernova (SN) 
feedback then naturally creates over-pressurized regions 
of hot gas that rapidly expand, flattening the potential 
at the center of galaxies, leading to an irreversible 
expansion in the DM orbits. Multiple bursts of star 
formation push more DM to larger radii, transforming 
an initi ally steep, cuspy DM d ensity profile into a cored 
p rofile (IGovernato et al. l l2012| ). 

iZolotov et all (|2012h are the first to examine how this 
impact of baryons on the DM structure affects the evo- 
lution of satellites around a Milky Way-massed galaxy. 
They demonstrate that SN feedback reduces the central 
DM densities of satellites with M* > 10 7 M Q before in- 
fall. After infall, the presence of a baryonic disk increases 
the mass loss rate via tidal stripping for all satellites 
(jZolotov et al.M2012t IPenarrubia et al.H2010D . This tidal 
effect is particularly strong for those satellites that en- 
ter with cored DM halos, further increasing the discrep- 
ancy in the central masses predicted by DM+baryon and 
DM-only simulations. In this Letter we examine the ef- 
fect that this evolution has on the observed properties of 
the satellite population in Milky Way-mass galaxies at 
z = 0. We then use our results to interpret the observed 
dynamics of the Milky Way and M31 dwarf Spheroidal 
(dSph) satellites. 

2. THE SIMULATIONS 

The satellite sam ple that we use here is drawn from 
IZolotov et al.1 ()2012h . and we refer the reader to that pa- 
per for full details of the simulations and satellite se- 
lection. Briefly, two halos with virial masses of 7 x 
10 n M Q and 8 x 10 n M Q at z = were run with 
and without baryons . The DM-o nly simulations were 
run with PKDGRAV (jStadell 120011 ). while the simula- 
tions with baryons were run with pkdgrav's N-Body 
+ Smoothe d Particle Hydrodyn amics (SPH) successor, 
Gasoline (jWadslev et al.ll2()() Hi. 

In the SPH runs, gas can cool via primordial and 
metal lines following IShen et all (|2010l ) . A prescription 
for cre ation and destruct i on of H2 is implemented fol- 
lowing iChristensen et aLl (|2012f) , and star formation is 
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TABLE 1 
Satellite Properties 



Satellite ID 
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Note. — Column (1) total stellar mass in the satellite. Column 
(2) corresponding V-band magnitude, based on the age and mctallic- 
ity of st ar particles using the Starb urst99 stellar population synthesis 
models (fVazqucz &: Lcithcrcr 2005). Column (3) mass in HI gas in the 
satellite. Column (4) maximum circular velocity, V max . 

a This mass of HI would remain undetected by current HI surveys. 

tied directly to the presence of H2 gas. A uniform UV 
background turns on at z = 9, mimicking reionization 
(|Haardt fc Madaul [20011 ). Supernovae deposit 10 51 ergs 
of thermal energy into the surro unding gas following the 
" blastwave" scheme d escribed in lStinson et al.l ()2006f ). 

IZolotov et al.l (|2012|) show that the z = satellite pop- 
ulations of both galaxies are in good agreement with the 
Milky Way (MW) and M31 satellite luminosity functions, 
from —8 < My < —15. Properties of the satellites in 
this luminosity range are listed in Table 1, including a 
few that are still gas-rich at z = 0. This range excludes 
from each simulation one brighter, Magellanic-like satel- 
lite. All but one of the gas-free SPH satellites studied 
here have been matched to a satellite counterpart in the 
DM-only run. In the remainder of the paper, we report 
kinematic values at 1 kpc for the simulations in order 
to both ensure convergence of the rotation curve, and as 
a useful comparison to the observations that have kine- 
matic data measured at < 1 kpc. 

2.1. Baryonic Evolution with Mass 

With the exception of three of the lowest mass ha- 
los, all of the satellite progenitors retain gas until in- 
fall, allowing their star formation histories (SFHs) to ex- 
tend aHe^s^jxnth^Jrrfah (hi^^ 

see IGrebel fc Gallaeherl l200l IDellenbusch eFaTI [20081 
iWeisz et al.l 1201 lh . However, the nature of the star for- 
mation varies with halo mass. Heating from the uniform 
UV background, combined with early star formation and 
SN feedback, removes a substantial amount of gas from 
halos below 1O 9 M , leaving them gas poor. These halos 
lose 10 — 90 times more mass in gas than stars formed 
by the time of their accretion. The remaining gas has 
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Fig. 1. — The total v c at 1 kpc lor the SPH satellites as a function 
of their \Z-band magnitude, at z = 0. The empty red squares and 
empty blue diamonds arc v c at the half light radii for MW and 
M31 dSphs, respectively. 

low surface densities and is inefficient at producing stars. 
The low level of star formation that can occur prior to 
infall is relatively constant rather than stochastic. 

Halos above 10 9 M Q , corresponding to halos with M* > 
10 7 M Q at infall, arc able to retain more of their gas for 
longer, allowing them to have extended SFHs. Unlike 
the lower mass satellites, the SFHs of massive satellites 
tend to be episodic and bursty. With their deeper po- 
tential wells, the more massive halos achieve an initially 
higher star formation rate (SFR) than the lower mass ha- 
los as gas cools in the central galaxy. The subsequent SN 
feedback from the burst heats the surrounding gas, shut- 
ting off star formation for a period of time until the gas 
can again cool and continue with another burst. These 
bursty SFHs lead t o DM core creation prior to infall 
(|Zolotov et al.|[2012l) . 

Many subhalos lose most of their gas after infall and 
are gas-free by z = (see Table 1)0 While some subhalos 
lose their gas nearly instantly at infall, some are capa- 
ble of retaining their gas for an extended time, and even 
having a low level of star formation. However, the SFRs 
become strongly suppressed after infall, and no subhalo 
continues to undergo the bursty star formation that con- 
tributes to DM core creation. Instead, the mass within 
1 kpc is substantially reduced when gas is stripped after 
infall. 

3. INTERPRETING OBSERVED DWARF SPHEROIDALS 

Below we discuss the observational properties of the 
Milky Way and M31 dSphs in light of results that show 
that baryonic effects can reduce the central DM masses 
of luminous satellites that infall to a ~ L* galaxy. 

3 Satellites that retain gas until 2 = m ay be artificially gas -rich, 
due to inefficient stripping of gas in SPH ifAgcrtz et al. 200^1. 
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3.1. Satellite Kinematics at Redshift 

Figure [T] shows the total v c at 1 kpc of the gas-free 
SPH satellites (dSph analogs, black circles) as a function 
of My. Plotted with the simulation data are observa- 
tional results for d Sphs with My brigh ter than —8 in 
both the MW from [Walker et alj (|2009[ ) and M31 from 
iTollerud et al.l (1201 21) and s upplemented by the compi- 
lation of lMcConnachid (|2012] ) for And II and And VI. For 

the observational data, v c (r 1 / 2 ) = Via 2 , where a is the 
line of sight stellar velocity dispersion, and n/2 is the 
half-light r a dius o f the dSph. iWalker et afl (|2009( ) and 
IWolf et all ()2010f ) have shown that assumptions about 
isotropy are minimized at r 1 / 2 , making the masses deter- 
mined at ri/2 the most robust. The half-light radii are < 
1 kpc for all of the observed dSph. Our measurements for 
the simulated gas-free satellites are upper limits to the 
values derived at typical half light radii, as the rotation 
curves of these satellites either peak or continue to rise 
at 1 kpc. 

Figure [T] shows that all of the v c values for the simu- 
lated gas-free satellites are less than 22 km/s, consistent 
with v c results for MW and M31 dSphs. In contrast, 
of the 12 matched DM-only counterparts to the SPH 
satellites, nine have 20 < v c (lkpc) < 34 km/s, grossly 
inconsistent with the MW and M31 dSph data. Previ- 
ous studies using DM-only simulations have found that 
subhalos simulated in a CDM-context are significantly 
more dense than observations of the MW dSphs suggest 
(iBovlan-Kolchin et alJ^OTTl [2011 IWolf fc Bullockl l2012l; 



lHavashi fc Chiball20lj . IZolotov et al.1 (|2012D . however, 
have shown that the central 1 kpc of dwarf satellites 
can be strongly affected by both SN feedback and tidal 
stripping, reducing their central DM densities. When 
these baryonic effects are modeled self-consistently, as 
they are here, the kinematics of simulated satellites are 
well-matched to the observed kinematics of dSph satel- 
lites. 

Figure [2] shows the difference in the DM contribution 
to v c at 1 kpc between the SPH runs used here and the 
identical DM-only runs that neglect baryonic effects, as 
a function of My of the SPH satellites at z = 0. It is im- 
portant to note that these results include both the effects 
of SN feedback on the DM distribution of the satellites 
prior to infall, and tidal stripping after infall. In order to 
quantify the change in DM mass that this combination 
produce in the central regions of satellites compared to 
the DM-only case, we fit a linear regression to the data 
points in Figure [2] The resulting fit is shown by the dot- 
ted line, and yields A(v C: ikp C ) = —10.47 — 1.35 x My. 
The y-axis range of Figure [2] excludes from the plot one 
outlier with A(w c ) ~25 km/s|3 

Figure [5] demonstrates the amount by which DM-only 
simulations over-predict the central DM masses of satel- 
lite galaxies, and hence why such simulations find that 
their satellites are inconsistent with the kinematics of the 
MW dSphs. An explicit example of the change quanti- 
fied in Figure [2] is demonstrated in Figure [3] The total 
rotation curve (DM and baryonic mass) for an SPH simu- 
lated subhalo is shown in Figure[3J which closely matches 

4 This outlier recently passed directly through the disk of its 
parent galaxy, which doesn't exist in the DM-only run, making 
A (i; c ) for this galaxy larger than in other satellites. 
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Fig. 2. — The difference in v c of the DM rotation curves at 1 kpc 
between the SPH and DM-only simulations, as a function of the 
V-band magnitude in the SPH satellites, for all satellites listed in 
Table 1. The dotted line is a linear fit to the points. 




Fig. 3. — The 2 = rotation curves of a simulated satellite and its 
DM-only c ounterpart. The v c f or Fornax is over-plotted, based on 
the data in I Walker et al. (2009). The combination of SN feedback 
(before infall) and tidal stripping (after infall) substantially lower 
the v c of the SPH satellite by z = 0, and is in good agreement with 
the observed v c of Fornax. 

the derived v c at the half light radius of Fornax, along 
with the rotation curve for the same subhalo in a DM- 
only run without baryons. The DM-only satellite clearly 
over-predicts the central mass of this satellite. For almost 
all of our satellites, the DM-only runs produce satellites 
with 2-4 times more mass in the central 1 kpc than their 
SPH counterparts. 

It is common when making predictions for the prop- 
erties of subhalos within the CDM model to associate 
the most massive DM-only subhalos at infall with the 
most luminous satellites at z = (e.g., IKoposov et al.1 
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[2001 iBovlan-Kolchin et al.l f20TH [2012h . i.e., an abun- 
dance matching technique. Indeed, the most massive 
subhalos in the DM-only runs in this work are matched 
to the most luminou s satellites in the SPH runs at infall 
(jZolotov et all 120121 ). However, SPH satellites experi- 
ence evolution that DM-only runs do not account for. It 
is therefore incorrect to assume that the central masses 
of DM-only subhalos should predict the observed masses 
in the inner regions of luminous subhalos. For example, 
the DM-only satellite shown in Figure [3] is a massive sub- 
halo that would be assumed to host a luminous dSph at 
z = 0. However, this assignment would over-predict the 
central mass of this satellite by factors of a few. 

3.2. The Scatter in v c 

A long-standing puzzle in the observational data is that 
some of the bright MW dSphs with My ~ —13 have 
similar v c values to some of the much fainter galaxies with 
My ~ — 8. In the field, galaxies display a clear trend of 
incre asing stellar mass with halo mass (e.g. JMoster et ah! 
|2012| ). yet satellites that span two orders of magnitude 
in luminosity seem to show no trend in their central v c 
with mass. The comparable velocities of Draco (My = 
—8.8), and Fornax (My — —13.4), for example, have led 
IPenarrubia et al.l (|20P8|) to conclude that Draco formed 
in a halo 5 times more massive than Fornax, despite being 
70 times fainter. 

The simulated SPH satellites shown in Figure Q] fol- 
low a similar trend to the observed MW and M31 dSphs; 
at z — simulated satellites with My = — 9 can have 
comparable kinematics to satellites with My = — 13. At 
infall, however, SPH satellites follow an increasing trend 
of stellar mass with halo mass, with M* oc M^ r . The 
strong trend seen at infall (and in the field) is erased in 
the subsequent tidal evolution of these satellites, result- 
ing in the scatter observed at z = 0. 

Since both SPH and DM-only satellites undergo tidal 
stripping, one might expect that SPH and DM-only satel- 
lites would display a similar scatter in v c at z = 0. How- 
ever, when DM-only satellites are matched to the stellar 
mass of their SPH counterparts, DM-only satellites still 
show a trend of increasing v c at 1 kpc with luminos- 
ity at z = 0, unlike satellites simulated with baryons. 
There are two effects that lead to this divergence be- 
tween SPH and DM-only satellites. First, the presence 
of a baryonic disk in the SPH runs (the DM-only runs 
don't have a disk) enhances the rate of mass lost at each 
pericentric passage for all SPH satellites. This effect is in- 
creasingly significant with increasing ti me after infall and 
for satellites on more ecc entric orbits (IPenarrubia et all 
120101 : IZolotov et"aLll2012[) . For example, our sample con- 
tains three low luminosity satellites (My > —11) that 
are on highly eccentric orbits. These SPH satellites un- 
dergo a greater reduction in their mass and central v c 
than satellites on circular orbits. This results in the large 
spread of v c values for simulated SPH satellites fainter 
than My ~ —11, similar to the observational data in 
Figure [TJ Second, the influence of a baryonic disk is 
especially strong for satellites with shallow DM density 
profiles, like the SPH satellites in our sample brighter 
than My ~ —12. This effect exacerbates the amount 
of central mass lost in the most luminous subhalos, and 
contributes to lowering the central v c of these satellites to 
under 25 km/s. These combined effects lead to satellites 



across the luminosity range that appear to have similar 
masses at z = 0. Without the effect of the disk, and with 
their steep density profiles, DM-only satellites experience 
overall less mass loss for a given orbit, and therefore tidal 
effects do not introduce a large scatter in their v c . 

4. CONCLUSIONS & DISCUSSION 

We have demonstrated that simulations that account 
for the effects of SN feedback and enhanced tidal strip- 
ping on satellites result in a satellite population whose 
kinematic properties match the observed properties of 
the Milky Way and M31 satellites. Our findings are 
in sharp contrast to studies using DM-only simulations, 
which over-predict the central masses of satellites in com- 
parison to observations. 

Our results predict: 

1. All dSph satellites brighter than My ~ — 12 should 
have central DM density slopes that are shallow (or 
cored) out to ~1 kpc, and total masses several fac- 
tors lower than predicted by CDM DM-only simu- 
lations. 

2. If a satellite fainter than My ~ —12 underwent 
substantial stripping (~90% of it's total mass, 
IPenarrubia et alj 120081 ) . it could have started off 
as a more luminous dSph with a DM core. 

3. The presence of cored density profiles at radii under 
~ 500 pc in satellites fainter than My ~ —12 is not 
ruled out by our simulations, as we do not resolve 
this region. 

4. Regardless of the density profile, the tidal effects 
of the disk can dramatically lower the central DM 
mass of any satellite, depending on infall time and 
orbital eccentricity. 

Most attempts to measure the central density slopes of 
dSph s have relied on usin g the spherical Jeans equations 
(e.g.. lGilmore et al.ll2007l) . but the mass and anisotropy 
of the stellar orbits are degenerate in the Jeans model, 
maki ng the results highl y dependent on adopted assump- 
tions ([Evans et al.ll2009l) Rec ent works (|Wolf fe Bullockl 
I2012t iHavashi fe Chibal 120121 ) have attempted to over- 
come this issue by searching for maximum likelihoods 
in parameter space, but still adopt Jeans modeling. 
Schwarzchild modeling avoids some of the assumptions 
inherent in Jeans modeling, and has been applied to For- 
nax and Sculptor, with the conclusion that both have 
cored DM density profiles (|Jardel fc Gebhardtl 1201 2t 
IBreddels etaLl [20121) . 

Alternatively, studying dSphs that have multiple stel- 
lar populations that span varying radii avoids any as- 
sumption about a mass model or isotropy, allowing a di- 
rect fit to two independent derivations of the mass. The 
central density slopes of Fornax a nd Sculptor have been 
determined with this method dWalker fc Pefiarrubia! 
20111: lAmorisco fc Evansl I2012t IBattaglia et alj 120081: 
Agnello fc Evansll2012D . with all studies concluding that 
both Fornax (with slope measured out to ~1 kpc) 
and Sculptor (measured interior to 500 pc) favor cored 
density profiles, and exclude cuspy density profiles at 
high significance. It has also been argued that the 
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fact that Fornax's globular cluster orbits have not de- 
cayed by dynamical friction requires a cored DM profile 
(|Goerdt et al.ll205aiCole et alJl2012l) . 

Higher resolution simulations will be necessary to pre- 
dict if small (on scales < 500 pc) DM cores exist in lower 
luminosity satellit es, though the observationa l evidence 
(e.g., in Sculptor, IWalker k, Peharrubiall201lj) suggests 
that they exist. We note that the substantial reduction 
in the overall normalization of the central DM densities 
of satellites due to the presence of a disk will lead to a re- 
duced DM annihilation signal. This suggests that deeper 
searches for DM annihilation could help to constrain the 
central masses of these objects. However, it is already 
clear that the overall concentrations of the MW dSph 
population are lower than predicted by CDM DM-only 
simulations (iBovlan-Kolchin et af1l2012t iWolf fc Bullock! 
120121: lHavashi fc Chiball2012t ). The lowered densities of 
our simulated satellites are consistent with these obser- 
vational results. 
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